Adsorbate-limited conductivity of graphene 
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We present a theory of electronic transport in graphene in the presence of randomly placed 
adsorbates. Our analysis predicts a marked asymmetry of the conductivity about the Dirac point, as 
well as a negative weak-localization magnetoresistivity. In the region of strong scattering, quantum 
corrections drive the system further towards insulating behavior. These results explain key features 
of recent experiments, and are validated by numerical transport computations. 
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Graphene (the two-dimensional allotrope of carbon pQ) 
offers a chemically stable platform to host various chem- 
ical adsorbates [2J |3J |U [5J [7] . Their presence strongly 
affects electronic transport, making graphene based de- 
vices suitable for chemical sensoring. A key observation 
in experiments on chemically functionalizcd samples is 
a marked asymmetry of the conductivity as a function 
of a back gate voltage, which is used to steer the sys- 
tem across the charge neutrality point (the Dirac point, 
which separates the valence band from the conduction 
band). The conductivity becomes symmetric only when 
the sample is annealed by a strong bias current which 
dislodges the adsorbates. 

In this paper we provide a theory, supported by numer- 
ical simulations, that explains this experimental feature 
for covalently bonded (chemisorbed) adsorbates [51I51HU]. 
Our model is based on the tight-binding description of 
electrons in graphene [TT]. Chemisorbed molecules are 
incorporated into this description as laterally attached 
additional sites, where the onsite and coupling energies 
are extracted from the band structure of a graphene sheet 
with regularly placed adsorbates. Each type of adsorbate 
introduces a characteristic local energy dependent scat- 
tering potential in the graphene, and suppresses the con- 
ductivity on one side of the Dirac point, with only a weak 
effect on the other. For example, adsorbed H + yields al- 
most insulating behavior in n-type graphene, while in 
p-type structures the conductivity is close to that of 
clean material. For OH~ the role of the bands is re- 
versed. When combined with scattering from a random 
Coulomb potential [HJ H3], the resulting conductivity 
traces (shown in Fig. [l| are consistent with the findings 
in experiments [Sj [5]. For a small adsorbate concentra- 
tion, these conclusions can be drawn from kinetic theory. 
For larger concentrations, we implement a recently pro- 
posed renormalization group (RG) analysis [TU [TS] to 
account for systematic quantum corrections to the con- 
ductivity resulting from multi-adsorbate scattering. Our 
predictions are in good quantitative agreement with the 
results of the numerical transport computations. 

In a graphene sheet, the carbon atoms are held to- 
gether via sp 2 -hybridized covalent bonds, while the elec- 
tronic transport takes place by hopping along ir orbitals 




FIG. 1: (Color online) Conductivity of graphene in the pres- 
ence of adsorbates (concentration n;) and Coulomb impurities 
(concentration m) as a function of charge carrier concentra- 
tion n e (carrier density n e /A c with A c = 3v / 3a 2 /4 the area 
per carbon atom), for various ratios rn/m = (J3/2iv)x where 
the dimensionless parameter j3 ~ 1 characterizes the scatter- 
ing strength of the Coulomb impurities [121 1131 125] . (a) H + , 
(b) OH - , with parameters from first-principle computations 
(the insets illustrate sample segments). The results are based 
on Eq. Q. 

which can participate in covalent bonding with adsor- 
bates. The electrons in the 7r-band of graphene with ad- 
ditional adsorbed atoms can be described using a tight- 
binding Hamiltonian 

ft = -7 X! ° ] i Cm + ^2 n ri, (la) 

U n = e^ n d n + j t (clj n + c dj,). (lb) 
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The first term of Eq. (la I corresponds to the Hamilto- 



nian of clean graphene, where c/ are annihilation opera- 
tors on sites of the honeycomb lattice, and the nearest- 
neighbor coupling strength 7 = 2hvu/3a (with bond 
length a = 1.42 A) determines the Fermi velocity Vd at 
the Dirac point. Each adsorbed molecule contributes a 
term of the form (jTbj) , where d n is the annihilation oper- 
ator on the adsorbate site and a n is the host position on 
the honeycomb lattice. The adsorbate density is param- 
eterized by rii/A c , where n$ is the adsorbate-to-carbon 
ratio, and A c — 3\/3a 2 /4 is the area per carbon atom in 
graphene. 

The model Hamiltonian (JTJ can be justified by first- 
principle calculations. Here we consider the adsorbates 
H + and OH~, chosen because of the presence of ambi- 
ent water in many experiments, and which also can be 
selectively driven towards the graphene via electric fields 
[3]. Using density- functional theory (DFT) [16], we find 
that H+ is described by £; = 0.66 7, 7^ = 2.27, while 
for OH - £j = —2.9 7, 7i = 2.37. These energies are 
comparable to the graphene hopping energy 7 = 2.6 eV, 
which necessitates the nonperturbative approaches em- 
ployed in this paper. Our DFT calculations also confirm 
that the adsorbates form covalent bonds to well-defined 
host sites. In the fully relaxed configuration, both adsor- 
bates are aligned in the vertical direction of the graphene 
sheet (for illustration, see insets of Fig. [TJ. Other energy 
scales (such as shifts of the graphene onsite energies and 
next-to-nearcst neighbor couplings) are small and can be 
safely neglected [17j [18] . 

The main building block of our analytical considera- 
tions is the derivation of the scattering amplitude due to 
individual adsorbates, which enters the collision term of 
kinetic theory and also features as input into the renor- 
malization group analysis that captures quantum correc- 
tions due to multiple scattering. 

In the first step we self-consistently eliminate the ad- 
sorbate sites from the Hamiltonian (JT|) via a decimation 
procedure. The electron wavefunction can be written as 

l*> = EiV»l 10 + E„^n |n) ad = m e where the 
amplitudes ipi and states \l) refer to the carbon sites, 
while 4> n and \n) ad refer to the adsorbate sites. We 
now project the Schrodinger equation onto an adsorbate 
site, (n| d (e — TL) Y$>) = and find that the amplitude 
<fin = V J Q„7i/( e — e i) on the adsorbate is related to the 
amplitude of its host carbon site. The amplitudes <p n 
hence can be eliminated, which results in the reduced 
Hamiltonian 
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n = -7 V c) Cm + Yvd c Q „, u = ^^, (2) 



where the energy-dependent effective potential V corre- 
sponds to the self-energy which an adsorbate induces for 
electrons in graphene. 

While the effective potential V displays a distinct res- 
onant energy dependence, the analysis of the resulting 



conductivity properties requires combining this with the 
specific energetics of the graphene sample, and in partic- 
ular, with the existence of the conical point at which the 
density of states drops to zero. We therefore now turn to 
the analysis for the Green's function Q = {e — H + i0 + )~ 1 
of the system. 

For a single, well isolated adsorbate, the Green's func- 
tion can be obtained exactly over the entire energy range 
by utilizing the T-matrix representation 

G = G + G TG , T = (l-Vclc a Go)- x Vclc a , (3) 

where Q$ is the Green's function of the clean graphitic 
system and T characterizes the scattering strength [T!5] . 
Expanding T in powers of VGq, and utilizing the trans- 
lational and crystalline symmetries of graphene (so that 
(n\Q Q \n) = g for all n), the resulting series for Q can be 
resummed yielding 



T = t (e)c< a c a , t (e) 
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i! 90(e)' 



(4) 



In the latter expression, to(e) describes resonant scat- 
tering of electrons in graphene from an adsorbate level 
renormalized by hybridization with states in the 7r-band. 
This hybridization always shifts the effective resonance 
level towards the Dirac point, which can be understood 
as a consequence of level repulsion which pushes addi- 
tional states towards the region with the lowest density 
of states. The real part of 'yfgo(^) [with Rego(e) = 
-Rego(-e) ~ (^ c /27rfi 2 u|))eln(|e|/A), and A a high- 
energy cut-off] gives a formal description of this sys- 
tematic shift, while the imaginary part [with Irng — 
—ttvq{e), and Uq(e) the density of states per carbon atom] 
indicates a decrease of the resonance width near the Dirac 
point. 

Ignoring (for the moment) the effects of phase-coherent 
multiple scattering off the adsorbates, the conductivity of 
a disordered sample with a finite adsorbate concentration 
rii can now be obtained in kinetic transport theory. Due 
to Cq rotational symmetry of graphene, the conductivity 
is isotropic. Using the scattering rate obtained from Q, 



"i|*o(£k)| 2 ^o(^k), 



(5) 



where £k is the graphene dispersion relation, we find the 
conductivity of graphene in the presence of the adsor- 
bates, 



9 s e 



h J A c m\to(E F )\ 2 (v^ 1 ), 



(6) 



where 



denotes integration along the Fermi line, 



Vk = |Vk£k|/7i, and g s = 2 accounts for spin degeneracy. 

Motivated by recent experimental findings for selec- 
tively adsorbed H + and OH - [3], we show in Fig. [2] the 
corresponding dependence of the conductivity in Eq. ^ 
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FIG. 2: (Color online) Conductivity of graphene in the pres- 
ence of adsorbates of variable concentration rij. (a) H + , (b) 
OH - . The dashed thick curve is the prediction m of kinetic 
transport theory with energy-dependent T-matrix The 
solid curves in the left panels show the expected quantum 
corrections based on the renormalization group analysis [Eqs. 
1 10 1]. In the right panels, the solid curves show the results of 
numerical transport computations. The inset in the left panel 



of (b) shows the RG flow ( 10 1 of the dimensionless scattering 
parameters for OH~ at e = —0.157 and rii — 0.01. 



on the Fermi energy (thick dashed curve). A key fea- 
ture is the marked asymmetry of the conductivity about 
the Dirac point caused by the energy dependence of the 
T matrix For H + , over a range of energies in the 
conduction band the conductivity is small, while it rises 
linearly as one moves into the valence band, or far into 
the conduction band. For OH - the role of the bands is 
reversed. 

In actual devices, adsorbate scattering is supplemented 
by scattering off localized charges. These contribute to 
the scattering rate a term t^ 1 = ni (3~/ 2 / (h\e\) , where n; 
is the number of localized charges per carbon atom, and 
j3 ~ 1 is a dimensionless numerical factor [T^l Q2]. For 
energies in the linear part of the clean graphene disper- 
sion relation, the resulting conductivity can be written 
in the form 



(x\t (n e )/-y\ 2 + n' e 



where x — 
amount of 



(2n/f3)(n t / ni ) 
the two types 



characterizes 
of disorder, 



) , (7) 

the relative 
and n e — 



f/ s £p/(2v / 37T7 2 ) is the number of charge carriers per car- 
bon atom. Figure [l] shows the charge carrier dependence 
of the conductivity for various values of x. For x = the 
conductivity shows the symmetric linear charge-carrier 
dependence characteristic for charged impurity scatter- 
ing. For an increasing adsorbate concentration, the 



conductivity develops the marked asymmetry discussed 
above, while the Coulomb scattering only dominates very 
close to the Dirac point, where it results in an additional 
dip. 

Different types of disorder can also be discriminated 
via mesoscopic corrections to the conductivity, which 
originate in phase-coherent multiple scattering. In gen- 
eral, disorder in graphene can be characterized in terms 
of five dimensionless parameters T = {a , P± , Z , 7j_ , 7z } , 
which classify the breaking of the symmetries of the hon- 
eycomb lattice [20 . Three of these parameters describe 
intravalley scattering preserving the Cq v symmetry (ao), 
the Cs v symmetry (j z ), or no point symmetry (7.1), 
respectively. Analogously, intervalley scattering is de- 
scribed by the parameters (3 Z and /3j_ . In terms of these 
parameters, the conductivity of graphene in kinetic the- 
ory takes the form 



Q s e f ao . a . 3 



+ o7; 



(8) 



The scattering potential for adsorbates is such that 



#l/2 



A c rn\to(e F )\ 2 

27T(v k ) £F /(v^) ( 



z = 1x = 0. 



This puts chemically functionalized graphene into the so- 
called orthogonal symmetry class, for which one expects 
a negative weak-localization magnetoresistance [5D], as 
well as strong Anderson localization when the adsorbate 
concentration increases [14] . 

These corrections to kinetic theory can be studied sys- 
tematically via a renormalization group analysis [21 , 
which provides effective scattering parameters T that re- 
place the bare values in Eq. ([8]). For disorder represen- 
tative of the symmetries in graphene, the flow equations 
were derived in Ref. |14j . An equivalent formulation in 
terms of the parameters T given above can be found in 
Ref. [IS], and takes the form 

d = 2a (a + f3±+j±+ f3 z + ~f z ) + f3 ± l3 z + 2jj_~f Zl 

fix = 4(«oA + /3±7-L+ Mz), 

f3 z = 2(a /3j_" (3 z a + f3±_l z + Mz), 

i± = 4a -f z +Pl + (3 2 z , 

i z = 2-f z (-a - /3l + P z + 7J- ~ 7*) + 2a 7-L + P±/3 Z , 
e = e(l + a Q + /3± + j± + p z + lz ), (10) 

where X = dX / d\sx{L / a) and L is a running length. 
This RG flow is integrated using the bare parameters T 
and e — 8S initial conditions and terminated when e 
reaches a high-energy cutoff e c . We here implement this 
procedure for the energy-dependent initial conditions (|9| . 
A typical solution of the flow equations as a function of 
hx{L/d) is shown in the inset in Fig. |2jb) (left panel). 
The dominance of /3± over 7 Z and 7j_ demonstrates that 
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chemically functionalized graphene stays in the orthogo- 
nal symmetry class when the adsorbate concentration is 
increased. As shown by the curves in the left panels of 
Fig. [2] over the whole energy range and for both types of 
adsorbates (H + and OH~) the renormalization leads to 
a suppression of the conductivity and drives the system 
towards insulating behavior. 

In order to verify these expectations we carried out nu- 
merical transport computations, which are directly based 
on the model Hamiltonian (|T|) . The conductivity a is ob- 
tained by finite-size scaling of the Landauer conductance 
G = aW/L of graphene ribbons with respect to width (in 
the range 50 < W/V3a < 100) and length (in an adap- 
tively chosen range that avoids the onset of Anderson 
localization [22]). For each geometry the conductance 
is computed in an efficient recursive Green's functions 
algorithm [23], and averaged over 10 4 disorder realiza- 
tions. The resulting dependence of the conductivity on 
the Fermi energy (right panels of Fig. [2]) is in good agree- 
ment with the expectations based on the RG analysis 
(left panels of Fig. [2J). In particular, the numerical re- 
sults clearly confirm the suppression of the conductivity 
due to quantum corrections. 

In summary, we have presented a theory of electronic 
transport for graphene in the presence of chemisorbed 
molecules. We find that each type of adsorbate results 
in a characteristic Fermi-energy dependence of the con- 
ductivity, which is asymmetric about the charge neu- 
trality point and distinguishes p- and n-type transport. 
In a range of energies, the conductivity is strongly sup- 
pressed, which could be used to increase the on-off ratio 
in graphene-based field-effect transistors. These effects 
are further enhanced by quantum corrections driving the 
system towards the localized state at higher adsorbate 
concentrations. In experiments, the adsorbate concen- 
tration can be increased by variable deposition times, 
or by driving adsorbates to the sample using top gates 
OUSEL while the localization effects can be probed via 
magnetoresistivity experiments on graphene flakes at low 
temperatures [24] . 
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